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Abstract. We study non-geodesic corrections to the quasicircular motion of charged test par- 
ticles in the field of magnetized slowly rotating neutron stars. The gravitational field is approx- 
imated by the Lense-Thirring geometry, the magnetic field is of the standard dipole character. 
Using a fully-relativistic approach we determine influence of the electromagnetic interaction 
(both attractive and repulsive) on the quasicircular motion. We focus on the behaviour of the 
orbital and epicyclic frequencies of the motion. Components of the four-velocity of the orbiting 
charged test particles are obtained by numerical solution of equations of motion, the epicyclic 
frequencies are obtained by using the standard perturbative method. The role of the combined 
effect of the neutron star magnetic field and its rotation in the character of the orbital and 
epicyclic frequencies is discussed. 
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We assume the external gravitational field of slowly rotating neutron or strange stars 
approximated by the Lense-Thirring metric (Lense & Thirring 1918, Hartle & Sharp 1967, 
Hartle 1967), the magnetic dipole with the symmetry axis identical with rotation axis and 
infinitely conductive star interior implying force lines frozen into the star and dragged 
by its rotation (Konno & Kojima 2000). Solving the equation of motion for a charged 
test particles, we obtain orbital angular velocities uj± for corotating and counter-rotating 
equatorial circular orbits (Fig. [I]). Properties of the corotational and counter- rotational 
circular orbital motion are naturally different due to the frame dragging of the external 
spacetime. The asymmetry is further amplified by the existence of two components of the 
Lorentz force with different behaviour. The Coulombic part of the interaction is induced 
by star rotation only, and is therefore independent of the test particles orbital velocity, 
while the orientation and amplitude of the magnetic part of the interaction induced by 
the orbital motion are strongly dependent on the orbital velocity. 

Surprisingly, in the case of the counter-rotating orbits, starting from a critical specific 
charge q os , the character of the electromagnetic interaction enables the existence of the 
electrostatic radius r es (q), where properly charged particles are static relative to static 
observers at infinity. In such case, the gravitational attraction of the neutron star is 
compensated by the electric repulsion (Fig. [I]) . 

The existence of epicyclic frequency corresponds to the condition of the stability rela- 
tive to the appropriately oriented perturbations. Therefore, vanishing of the radial u r or 
vertical ujq epicyclic frequencies determines the location of the marginally stable circular 
orbits (Fig. [2j [3]) . We use analytic formulae for the radial and vertical epicyclic frequen- 
cies of a charged test particle in the presence of a general electromagnetic field derived 
and discussed by Aliev & Galtsov 1981, Aliev 2007, Bakala et al. 2010 k 2012. 
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Figure 1. Contour plot of the orbital fre- 
quency v — ui/2-n as a function of the spe- 
cific charge q and the radial coordinate. Top 
left: Corotating solution for spin a = 0.05. Top 
right: Corotating solution for spin a = 0.3. 
Bottom left: Counter-rotating solution for spin 
a = 0.05. Bottom right: Counter-rotating solu- 
tion for spin a = 0.3. 




Figure 2. Top left: Contour plot of the coro- 
tating radial epicyclic frequency v r ± = cu r ±/2n 
as a function of the specific charge q and the 
radial coordinate. Bottom left: Contour plot 
of the vertical epicyclic frequency ve = uug/2n 
as a function of the specific charge q and the 
radial coordinate. Right panels: Same as the 
left panels, but for counter-rotating orbits. 
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Figure 3. The region of stable circular or- 
bits filled up by the contour plot of the or- 
bital frequency v = Q./2-k. Constructed for 
a test neutron star with M = 1.5M© and 
fj, — 1.06x 10 -4 m 2 . Top left: Corotating or- 
bits for spin a = 0.05. Top right: Corotating 
orbits for spin a = 0.3. Bottom left: Coun- 
ter-rotating orbits for spin a = 0.05. Bottom 
right: Counter-rotating orbits for spin a = 0.3. 

Acknowledgments 

This work has been supported by 
the Czech grants GACR 202/09/0772 
and GACR P209/12/P740. The authors 
further acknowledge the internal stu- 
dent grants of the Silesian University in 
Opava (SGS/1/2011, SGS/2/2011) and the 
project CZ. 1.07/2.3.00/20.0071 "Synergy" 
in the frame of "Education for Competi- 
tiveness" Czech Operational Programme. 



References 



899 



Stuchh'k Z 
065012 



Aliev A N & Galtsov D V, 1981, GRG, 13. 

Aliev A N, 2007, |astro-ph/0612730"V 2 

Bakala P, Sramkova E, Stuchh'k Z and Torok 
G, 2010, CQG, 27, 045001 

Bakala P, Urbanec M, Sramkova E 
and Torok G, 2012, CQG, 29 ; 

Hartle J B, 1967, Ap.J., 150, 1005 

Hartle J B & Sharp D V, 1967, Ap.J., 147, 317 

Konno K & Kojima Y 2000, Progress of The- 
oretical Physics, 104/6, 1117 

Lense J und Thirring H, 1918, Phys. Zeitschr, 
19 



